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I. INTRODUCTION 

It has been recently pointed out ^ that, albeit integrals containing products of Gaussian 
unctions and of an arbitrary number of Hermite polynomials are available since long time 



2|, the same is not true 
has been put forward in 



or Laguerre polynomials. A first step aimed at filling this gap 
J|, and in this paper we give our contribution using a recently 
developed js] symbolic method, which has proved its usefulness in the study of integrals 
involving Bessel functions. 

Following this symbolic approach, the two-variable Laguerre polynomials 



. k „,n—k 



^ [n-k)\[k\Y 

can be written as 



Ln{x,y) = {y - cxY ipo (1.2) 
where the following umbral notation has been used [s], Q] 

^>o = ^- = f(^- ^'-'^ 
^Pq is referred to as the polynomial vacuum (the umbral operator c acting on the "state" </?o 

"creates" the Laguerre polynomials). 

Before closing this short introduction, we report a result that will be useful in what 

n 

follows. As it has been done in Ref. [5|, by using the generating function of the two- variable 
Hermite polynomials 

o° j-n n-2k k 

j:^H^(.,y) = e^'-"' iJ„(,,,) = „!5]^_^, (1.4) 

n=0 k=0 ^ ' 

it is possible to show that 

/„(..M.ffl=£d.(o. + .)V-"^'«^ = y^exp(g //„(f+^.|^) . (1.5) 

The use of the symbolic method along with integrals of the type (11. 5p will be the key note 
for the future development. 

II. A SURVEY OF INTEGRALS OF LAGUERRE POLYNOMIALS 

By taking into account the expression (II. 2p for the Laguerre polynomials, one can write 
Xn{u;a)= / dxLn{x,u)e ""^ = /„,(-c, w; a, 0) v?o , (2.1) 



and, treating c as an ordinary constant, it's easy to show that 

^.(";«) = y^«f (2.2) 

where we have introduced the polynomials 

Ql:\x,y)=c''H^{x,yc^)ipo 

^n-2kyk 

" ^ (n-2A;)!A;!r(2A; + z/ + l) ' ^^'^^ 
whose properties will be discussed later in the paper. 

The form fll.2p generalizes to the case of associated two-variable Laguerre polynomials as 
follows 

Ll:\x, y) = + c^y-cxT (2.4) 

and, proceeding as above, we get 

xr)(.;«)^£d.Lr(^.")e--^;^^<^i±i^e'(".J^) . (2.5) 

These results and the identity^ 

Hr,{x + a,y) = Y, Cf) Hn^kix, y) (2.6) 
allow us to obtain the following formula 

/oo 
dx L^^\x + a,u) e"""' 

= \la n\ ^ U ) ^-"^ ^"-^ 4^J ■ ^ ^ ^ 

The multi-index extension of the Hermite polynomials can be exploited to compute in- 
;egrals involving products of powers of binomials and gaussians in a compact form. In Ref. 
6| it has been shown that 



r dx (a X + hr ifx + gre ""^^ = . ( g, ^] , 
j_aa \ Oi \ 4: a 2aJ 



(2i 



^ It is essentially a Taylor series expansion, and can be proved by noting that 

oo , fx \k 

U^(x + a, J/) = e"^"^ y)^Y. ^"(^' y) 

k=0 

and 
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where 



min(m,n) 

Hm,n{x,y;w,z\T) = m\n\ ^ 



k=0 



{m — k) \ {n ~ k) \ k\ 



Hm-k{x,y) Hn-kiw,z) (2.9) 



are the two-index Hermite polynomials, whose properties are discussed in detail in Refs 



and, more recently, in 



.fl 



J] . The use of the same argument as before yields the following result 



n r(m + + T{n + u + l)_ 



m! 



Aa Aa 



1 
2a 



(2.10) 



where 



min(m,n) 

Qt%\^^y^w,z\T) =m\n\ 



k=0 



{m - k)\ {n - k)\ k\ 



Qtt^{x,y)Q^::^;:\w,z). (2.11) 



(The case of the product of two Laguerre polynomials is obtained putting /i = z/ = 0). The 
result for an arbitrary product of associated Laguerre polynomials can be easily derived from 
the previous identities. It involves multi-index Q-polynomials and it is not reported for the 
sake of brevity. 

The method is flexible enough to be easily extended to the calculation of family of integrals 
containing products of Laguerre and Hermite polynomials. Taking into account Eq. (9) of 



Ref. lol, one obtains 



dx (a X + br Hnif x + g,y) e""^'' = ^^^-.n ^■^9,y + 



Aa 



af 
2a 



(2.12) 



and, thus, by using expression (11.21) for the associated Laguerre polynomials, it's easy to 
show that 

^ r(n + i^ + ll 



r dx (x, ,) H^if x + g,z) e—' = .[^ + ^ + (y,^;g,z + f-I-) 

J_^ \ a nl ' \ Aa A a la) 

(2.13) 



where 



min(m.,n) ^ 

rS(x,,;..,z|r)=m!n! g _ fc),^ - fc)! fc! 
As a final example, let us consider now the following integral 

/oo 
dx Ln{x,y) Jo{2V^) e"""^ 
-oo 



Q^J:,t:\x,y)R,.k{w,z). (2.14) 



(2.15) 



where Jo{x) is the first-order cyhndrical Bessel function, that, according to Ref. can be 
cast in the form 

Jo(x)=exp|-d (^)'l Ao (2.16) 

with the umbral operator d acting on the polynomial vacuum Aq in the same way that c acts 
on Lpo (see Eq. (11. 3p ). This integral looks quite complicated, but using the expression (II. 2p 
for the Laguerre polynomials, and taking into account Eq. (II. 5p . we can write 



Br^iy, a) = a/^ ( 1/ + 1^ ) exp ^ ^ v^o Ao 



Since 



we get 



where 



X 



Bn{y,a) 



k=0 
TV 

a 



k\T{2k + iy + l] 



k=0 



n-2k 



An-2k{u, V)=}_^[^, 1 ^ ^ W, 



p=0 



2k + p 



Av 



(2.17) 



(2.18) 



(2.19) 



(2.20) 



The function VFq,(x|2) is the second-order Bessel- Wright function [81], that intervenes in the 
expression of the generating function of the Q-polynomials defined in Eq. fl2.3p . i.e. 



V ^ Ql^\x, y) = c° e^'+y'' '\o = e^' W^{ye\2) 



(2.21) 



n=0 



where the expression for the generating function of the Hermite polynomials (see Eq. (11. 4p ). 
and the identity (I2.18P have been used. This means that Q-polynomials belong to the Sheffer 
family and can be defined by the operational identity 



Q^:\x,y) = W^{y6i\2)x'^. 



(2.22) 



They cannot be considered orthogonal polynomials in a strict sense of the word, but pseudo- 
orthogonal according to the notion introduced in 9[. This aspect of the problem will be 
addressed in a forthcoming publication. 
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